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^ ■ Abstract 

We study cosmological evolutions of the generalized model of nonlinear massive gravity in which 
the graviton mass is given by a rolling scalar field and is varying along time. By performing 
dynamical analysis, we derive the critical points of this system and study their stabilities. These 
critical points can be classified into two categories depending on whether they are identical with 
the traditional ones obtained in General Relativity. We discuss the cosmological implication of 



relevant critical points. 
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I. INTRODUCTION 



Initiated by Fierz and Pauli (FP) [If, it has been questioned for long time whether 
the graviton is allowed to acquire a mass and leads to a consistent covariant modification 
of General Relativity. At quadratic order the FP mass term is the only ghost-free term 
describing a gravitational theory containing five degrees of freedom [2], but this theory can 
not recover linearized Einstein gravity in the limit of vanishing graviton mass, due to the 
existence of the van Dam-Veltman-Zakharov (vDVZ) discontinuity arising from the coupling 
between the longitude mode of the graviton and the trace of the energy momentum tensor 

3|, |4|. It was later noticed that this troublesome mode could be suppressed at macroscopic 
length scales due to the so-called Vainshtein mechanism jlj] . However, these nonlinear terms, 
which are responsible for the suppression of vDVZ discontinuity, lead inevitably to the 
existence of the Boulware-Deser (BD) ghost |6j] and therefore make the theory unstable 

7H1Q|. 

Recently, a family of nonlinear extension on the massive gravity theory was constructed 
by de Rham, Gabadadze and Tolley (dRGT) 
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121 ] . In this model the BD ghosts can be 



removed in the decoupling limit to all orders in p erturbation theory through a systematic 
construction of a covariant nonlinear action 
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16( | (see [l7j for a review). As a conse- 



quence, the theoretical and phenomenological advantages of the dRGT model led to a wide 



investigation in the 



are discussed 



40 



in [l8 



iterature. For example, cosmological implications of the dRGT model 



39| | ; black holes and spherically symmetric solutions were analyzed in 



50] ; and connections to bi- metric gravity models were studied in 5lM64|. 



Among these phenomenological studies, a generalized version of the dRGT model was 
constructed in Ref. that the graviton mass can be determined by a rolling scalar field. 
This mass-varying massive gravity (MVMG) model was argued to be free of the BD ghosts 
as well through examining the constraint system in the Hamiltonian formulation. Thus, it is 
interesting to investigate the cosmological implications of this model, especially its late time 
evolution. In the present work we perform a phase-space and stability analysis of this model, 
and investigate the possible cosmological dynamics in a systematical way. This method was 
widely applied in the literature and was proven to be very powerful particularly in the study 



of dark energy physics 



65 



476] (see also Refs. 73, |78|] for relevant reviews). In the present 



model the cosmological system shows a couple of critical points at late times, and they can 
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generally be classified into two categories depending on whether they are identical with the 
traditional ones obtained in General Relativity (GR). Although some of the critical points 
are able to be recovered in GR in the limit of vanishing graviton mass, their background 
dynamics are different since of the graviton potential. By performing stability analysis we 
find that the parameter space of the model is tightly constrained and the stable cosmological 
evolutions are quite trivial. Moreover, there exist some new critical points in the MVMG 
model which might be of theoretical interests, but the corresponding cosmological evolutions 
are basically ruled out by observations. 

The present paper is organized as follows. In Section II, we briefly review the MVMG 
model and its cosmological equations of motion. Then we perform a detailed phase-space 
and stability analysis of this cosmological system and summarize the results in Section III. 
Finally, we conclude with a discussion in Section IV. 



II. THE MVMG COSMOLOGY 



To begin with, we briefly review the MVMG model constructed in 35(. This model 
requires the graviton mass to a function of another scalar field by introducing a nontrivial 
potential term V(ip) and thus it is able to vary along cosmic evolution. General Relativity 
can be automatically recovered when this scalar field ip sits at V(ip) = 0. In order to better 
control dynamics of the scalar field, the model also includes an additional potential W(ijj). 
Therefore, the complete action can be expressed as 



S 



d %\ 



-fj 



M 2 i 
-^■R + V (VO (U 2 + a 3 U 3 + a A U 4 ) - -d^d^ - W (V) 



(1) 



The U21 U3 and U4 terms are the graviton potentials following the dRGT model, which are 
given by, 



with 



U, 



K>t = 8$- \ gwf AB d p <j) A d v <t> B . 



Ua 



(2) 



(3) 



Moreover, f^B is a fiducial metric, which is often chosen as Minkowski Jab — Vab- The four 
4> A (x) are Stiickelberg scalars introduced to restore general covariance under 



dx p dx c 
dx'v dx" 



■g pa {x) , <p A (x) -> <p A (x) 



x* ->■ x'" 



(4) 



For the case Jab = Vab, the Stiickelberg scalars form Lorentz 4- vectors in the internal space. 
By performing the Hamiltonian constraint of the system it is argued that this model is still 



free of BD ghosts [351 ] . 



In the following we take into account the regular matter component in the total action 
which minimally coupled to the gravitational system. Further, we consider the fiducial 
metric to be Minkowski 

Sab = Vab, (5) 

and assume that the dynamical and fiducial metrics are of diagonal forms simultaneously 
for simplicity. Specifically, we consider a flat Friedmann-Robertson- Walker (FRW) metric 

d 2 s = -N(r) 2 dr 2 + a(r)%dxW, (6) 

and for the Stiickelberg fields we choose the ansatz 

(f>° = b(r), (t> l = a x\ (7) 

where a is constant. 

Finally, the total action for the cosmological background reduces to 

S T = I d 4 x 



+ S„ 



where 



u 2 = 3a(a — a )(2Na — a N — ab'), (9) 



u 



f _ I „ „ \2 
3 



(a-a Y{ANa-a N -3ab'), (10) 



v% = (a-ao) 3 (N-l/), (11) 

and we have define ' = d/dr. Variations of the total action St with respect to N and a 
lead to the two Friedmann equations 

3M 2 P H 2 = PMG + Pm, (12) 

-2MpH = p MG + p MG + pm+Pm , (13) 

respectively. In the above equations, p m and p m are the density and pressure for the matter 
component, respectively. The effective density and pressure are given by 

Pmg = \tf + Wty) + VW){X - l)h{a h X) + VW)(X - l)h{a h X), (14) 

Pmg = - Wty) - V(iP){X - l)f 3 (oci,X) - V(iP)(b - l)fi(oi,X), (15) 
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with 



h{a h X) = (3 - 2X) + a 3 (3 - X)(l - X) + a 4 (l - X) 2 , (16) 



/ 2 (ai, X) = (1 - X) + a 3 (l - X) 2 + ^(1 - X) 3 , (17) 
f 3 ( ai ,X) = (3-X) + a 3 (l-X), (18) 
and X = a /a. 

Moreover, variations of the total action with respect to b and ip give rise to the following 
two equations of motion, 

V{i))Hh{a h X) + V{i))f 2 {a h X) = , (19) 
dW dV 

^ + 3H^ + —- + —[(X- l)(f 3 ( ai , X) + /!(«,, X)) + 36/ a (ai, X)] = . (20) 

which will be frequently used in detailed calculations later. In addition, it is convenient to 
define an effective equation of state parameter for modified gravity terms as follows, 

Pug W + WW + ViMX - l)h{a h X) + V{^){X - l)h{a<, jQ , 01 , 

Wmg = = ^r~- = ; (21) 

Vmg ^-W(i>)-V(i;)(X-l)f 3 (a i ,X)-V(i;)(b-l)f 1 (a i ,X) 

and the total equation of state parameter is defined as: 

PMG + Pm 1 2H 

Vmg + Pm -J-n 
III. DYNAMICAL FRAMEWORK OF MVMG COSMOLOGY 



In this section we perform a detailed phase-space analysis of cosmic evo 



utions described 



65M75|| (see also 



by the MVMG model. Following the method extensively developed in 
for a recent analysis in the frame of generalized Galileon cosmology), we first transform the 
dynamical system into the autonomous form. 



A. Dynamics of the autonomous system 

In general, for a dynamical system one can suitably choose a group of auxiliary variables, 
and the corresponding equation of motion can be expressed as a group of first-order differ- 
ential equations, respectively. Namely, to illustrate the method of phase-space analysis, we 
consider the following two- variable dynamical system: 

i = ffay) , y = g(x,y) . (23) 

5 



The system is said to be autonomous if / and g do not contain explicit time-dependent terms. 
A point (x c ,y c ) is said to be a critical point of the autonomous system if (/, g)\f Xc , yc ) = 0. 
One can check whether the system approaches one of the critical points or not by performing 
the stability analysis around the fixed points. Specifically, one can introduce Sx and Sy as 
small perturbations and expand the differential equations ( 12 3 j) to the first order of 5x and 
Sy around the critical point, and then can derive out the following equations of motion, 




dx dy 
dg_ dg 
dx d y / (x=x c ,y=y c ) 




(24) 



As a consequence, the general solution for the evolution of linear perturbations can be 
written as 

5x = Cl e» lN + c 2 e" 2iV , (25) 



Sy = c 3 e» lN + c 4 e" 2iV , (26) 

where [L\ and /i 2 are the two eigenvalues of matrix in the left hand side of Eq. ( 124)) . and 
Ci, c 2 , c 3 and c 4 are constant coefficients. If < and fi 2 < 0, then the point is stable, 
which means the system could evolve to this fixed point eventually. The method can be 
extended to a system with many variables, a critical point is stable if the real parts of all 
the corresponding eigenvalues are negative. 

In the model we consider, there are five dimensionless variables 

Xp ~ V3M P H' Xf ~ V6M P H' XW ~ V3M P H' %V ~ V3M P H' %a ~ a ' 
Among them, x a and Ob p CcLll be determined by background equations of motion as will be 
analyzed in this subsection. Thus the system only involves three independent variables. 
Making use of these variables, one can reexpress the Friedmann equation as follows, 

1 = x p 2 + x/ + x 2 w + x v 2 {x a - l)(/i(0 + /sW) • (28) 

which now is a constraint equation. 

Then, we particularly choose an exponential potential 

A 



W = exp 



-ib 



(29) 



with A > 0. Moreover, we parameterize the form of b as a linear function of cosmic time, 
which is given by 



b = Bt with B > , (30) 

so that this autonomous system is analytically solvable. In addition, we also assume the 
matter fluid satisfies a barotropic equation of state p m = (7 — 1) p m , with 7 being a constant 
and < 7 < 2. From equation ( |T9l . one can get 

= (31) 

dip / 2 \X a ) 



and thus 

dV VHh (x a ) 



# i>h (X a ) 



(32) 



if ip ^ 0. Note that, if ^ = the mass term Vjjp) is fixed and then the model would reduce to 



the dRGT version which has been shown in 



19| that a flat FRW background is not allowed. 



Therefore, we will not consider this case in the present work. Moreover, using the auxiliary 
variables, one can transform the equations of motion to the following autonomous forms, 

f-Xp {^fX 2 + 2x1 + fl ^a) {X a ~ B) Xy) - ~7X p , (33) 



Hdt Xp 2~ PKI ~ P ' ~~ w ' ^~«/v~« 2 

+ x v f 1 (x a ) _ + + wh{Xa)) (34) 

2x^f 2 {x a ) 

B^lkt Xw = ~Y~^ lx2 p + 2x i + f^ Xa ^ x v^ Xa ~ B )) ~ ~y- x ^ x w , (35) 

iijt xv = ^ (7x ' + 2x ^ + f^ x v^ - b » - ' (36) 

Hlt Xa = ~ Xa ■ (37) 

We restrict our discussion of the existence and stability of critical points to the expanding 

universes with H > 0. The critical points correspond to those fixed points where j^x p = 

0' lhtt x ' l l ) = 0' ~Wdi Xw = 0> ~Wdi Xv = an d lhi Xa = 0> an< ^ there are self-similar solutions 
satisfying 

H 3 

= ~2^4 + 2x l + h( x a)x v (x a - B)). (38) 
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The equation ( 157)) simply suggests x a = 0. It can be shown that having the variable x a 
in the system always bring a eigenvalue of —1 and leave the other eigenvalues unchanged. 
Thus in the following, for simplicity we will take x a = 0, which dose not affect the final 
results. 

By defining 

f x = 3 + 3a 3 + a 4 ,/ 2 = y,/ 3 = 3 + a 3 , (39) 
and using x a = 0, the constrain equation reduces to 

1 = x p 2 + x/ + x 2 w - x v 2 {fi + / 3 ). (40) 
Further, we use the constrain equation to eliminate the variable x p and get 



Av/6 



2^(7(1 - x 2 w ) + (2 - 7 )zJ + cx v 2 ) - 2x^ + ^x 2 w -d^ = 0, (41) 

3x H /(7( 1 - ^w) + ( 2 ~~ l) x % + cx y) ~ ^VQx^x w = , (42) 

x y ( 7 (l - x 2 w ) + (2 - 7 )xJ + cx 2 v ) - x v = , (43) 
and the self-similar solutions reduce to 

3 

— = --( 7 (1 - x 2 w ) + (2 - 7 )xJ + cx 2 y ) , (44) 

where in order for convenience we have defined two dimensionless parameters c = (fi + 
/ 3 ) 7 - J B/ 1 andd=/ 1 (l-S) + / 3 . 

B. Phase-space analysis and results 

We summarize the fixed points of this autonomous system and their stability analysis in 
Table I and Table II, respectively. In the following we discuss these solutions case by case. 

From Table I, one may notice that the points are either of xy = or not. For those 
with xy = (which are the points (al), (a2), (b), and (c), respectively), the final state 
of the universe corresponds to that the background evolution is mainly determined by the 
exponential potential instead of the graviton mass. These solutions are able to be recovered 
in the limit of vanishing graviton mass and thus are consistent with the standard results 

8 



2 



Points that can be recovered in GR 



Lable x\y xy 

(al) 1 

(a2) -1 

76 V 1 "^ 



(b) -fc 



Points that can't be recovered in GR 
Lable x^ xw xy 



(d) y feg (j - 1) 

fell ../ d(7-1) / 1-7 



( e2 ) V c+d(7-2) °_ \/ c+d( 7 -2) 

TABLE I: The critical points in a spatially flat FRW universe with exponential potentials in the 
MVMG cosmology. 



Points that can be recovered in GR 



Lable 


Existence 


Stability 


Equation of state 


(al) 


all 7 and A 


unstable 


1 


(a2) 


all 7 and A 


unstable 


1 


(b) 


A 2 < 6 


stable if A 2 < min(37, 3) 


1 ^ 3 


(c) 


A 2 >3 7 


stable if 7 < 1 


-1 + 7 



Points that can't be recovered in GR 



Lable Existence Stability Equation of state 

stable if 7 > 1, A 2 > 3,E > 
stable if 7 > 1, A 2 > § (l + ^ry) 
unstable 

TABLE II: The properties of the critical points. The definitions of E is given in the text below. 

obtained in GR. A further stability analysis suggests the points (al) and (a2) are unstable. 
For the point (b), the solution is allowed if A 2 < 6 and stable against perturbations when 
A 2 < min(37,3). For the point (c), the solution exists when A 2 > 37 but is stable only if 



(d) >0,d>-j£ 
(<*) ^<0,^2j>0 

(«2) d <0'^fe L 2)>° 
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7 < 1. We list the eigenvalues of these solutions in detail as follows. 

Point (al) is a kinetic-dominated solution. The linearization of the system yields three 
eigenvalues 

3 l~3 
mi = -,m 2 = 6 - 37, m 3 = 3 - J -X . (45) 

The corresponding effective equation of state for the whole system approaches to w — 1 
which implies a stiff fluid dominant phase at late times of the universe. 

Point (a2) is a kinetic-dominated solution as well. The linearization of the system yields 
three eigenvalues 

3 l~3 
mi = -,m 2 = 6 - 37, m 3 = 3+ J -X , (46) 

and correspondingly, the destiny of the universe is the same as Point (al). 

Point (b) is a scalar-dominated solution since x^ 2 +xw 2 = 1, and only exists when A 2 < 6. 
The linearization of the system yields three eigenvalues 

mi = i(A 2 - 3) , m 2 = i(A 2 - 6) , m 3 = A 2 - 3 7 . (47) 

The corresponding total effective equation of state is w to t = — 1 + 4^ which is always less 
than unity. 

Point (c) is a solution depending on both the scalar field and the matter fluid, since 
x^ 2 + xw 2 = fi- The linearization of the system yields three eigenvalues 




The total effective equation of state is given by Wtot = — 1 + 7 in this solution. 

Point (d) is a solution depending on both the scalar field and the graviton mass since 
x % + x w ~ x v(h + ^3) = 1) the linearization of the system yields three eigenvalues 

mi = 3(1 - 7) , 

m 2 = -E+ ^3 (3 - A 2 ) xly + E 2 , 

m 3 = -E- ^3 (3-X 2 )x 2 w + Ei , (49) 
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, d('3-A 2 )(l-7) 

where E is a dimensionless parameter with E = | H — 2(0-^7) — ■ 

Point (el) is a solution determined by the kinetic term of the scalar field and the graviton 
mass since xyy = at late times. Its linearization yields three eigenvalues 




d ^- 1] (50) 



c + rf( 7 -2) ' 

Point (e2) is similar to Point (el) with its eigenvalues given by 



mi 
m 2 




It is clear that m 3 is positive, so the solution is unstable. 

It is worth noticing that the total effective equation of state of Points (al), (a2), (b), and 
(c) are the same as the massless case. This result shows that the effect of graviton mass does 
not contribute manifestly in the solutions corresponding to GR. Moreover, if we substitute 
the critical points (d), (el) and (e2) into the Eq. (I44j) . then we get H/H 2 = —3/2 which 
implies a matter domination at late times. This conclusion suggests that the appearance of 
a scalar field dependent graviton mass could strongly fix the background dynamics of the 
universe which manifestly conflicts with the observational fact of late time acceleration. 

We would like to point out that the existence of points (d), (el) and (e2) does not conflict 
with the fact that massive gravity theory does not allow flat FRW cosmologies, the fixed 
value of xy does not mean that the graviton mass is fixed. What actually happens is that 
the graviton mass and H both gradually approach as t — > 00 with the magnitude about 
1/t 2 and 1/t respectively. 

Finally, we can conclude that if a MVMG model is of cosmo logical interest, then its 
model parameters have to satisfy either A 2 < min(37, 3) (required by the stability of Point 
(b)) or 7 < 1 (required by the stability of Point (c)). Specifically, the solution of Point (b) 
corresponds to that the final evolution of the universe is determined by the scalar field and the 
effect of graviton mass is totally negligible, and thus this solution is quite trivial. Moreover, 
the solution of Point (c) corresponds to that the destiny of the universe is determined by the 

11 



combined effects of the scalar field and the matter fluid. However, the stability of Point (c) 
requires the barotropic equation of state of matter fluid to be less than unity, which implies 
the corresponding matter fluid has to violate the strong energy condition and thus obviously 
conflict with cosmological observations. 



IV. CONCLUSIONS 



In the present paper we have studied the dynamical behavior of the MVMG cosmology. 
This model, due to the varying of the graviton mass, possesses plentiful phenomenological 
properties and consequently has attracted many attentions in the literature. Namely, it 
was shown to be able to violate the null ener gy condition effectively and thus could be 



; further, it was also app 



ied into early 



applied to realize the phantom divide crossing 
universe and a class of bouncing and oscillating solutions were reconstructed [37| . However, 
while phenomenological studies of this model is still proceeding extensively, it is necessary 
to investigate the phase space of the model and examine the stability of the critical points 
existing in cosmological trajectories. Thus we performed a detailed dynamical analysis of 
the MVMG model with a specifically chosen potential for the cosmic scalar field, namely an 
exponential potential. Our result reveals that there are mainly two types of critical points 
in this model. One type of critical points correspond to the case that they can recover the 
standard results in GR if the graviton mass is chosen to be vanishing; the other type of 
critical points then is discontinuous with GR in the massless limit. We analyze these points 
respectively and find that there are only two critical points which might be stable against 
perturbations and both two belong to the first type. However, one of these two solutions 
is difficult to accommodate with cosmological observations since its stability requires the 
matter components in the universe to violate strong energy condition. Eventually, there is 
only one viable solution in this model but the final state of the universe is completely 
determined by the cosmic scalar field and there is no effect of modified gravity. As a 
consequence, the MVMG cosmology severely degenerate with the standard cosmology based 
on GR. 

We would like to point out that, although the graviton mass does not give rise to ob- 
servable effects on the background evolution at late times, it may still leave signatures on 
cosmological perturbations and thus affect the large scale structure formation. Moreover, in 



12 



our investigation, we focused on a particularly chosen potential for the cosmic scalar field. 
It would be interesting to generalize the case to a much more generic potential and verify if 
the phase space of the viable solutions could be enlarged. 

Note added: While this work was being finalized, we learned of a related work by 
G. Leon, E. N. Saridakis and J. Saavedra which will be appeared on arXiv 79j. Part of 
their content overlaps with ours, and their conclusions are similar as well while their focus 
is on a generalized structure of the phase space. 
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